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Abstract 

We study the dynamics of a charged tracer particle (TP) on a two-dimensional lattice all 
sites of which except one (a vacancy) are filled with identical neutral, hard-core particles. The 
particles move randomly by exchanging their positions with the vacancy, subject to the hard- 
core exclusion. In case when the charged TP experiences a bias due to external electric field 
E, (which favors its jumps in the preferential direction), we determine exactly the limiting 
probability distribution of the TP position in terms of appropriate scaling variables and the 
leading large- n (n being the discrete time) behavior of the TP mean displacement X n ; the 
latter is shown to obey an anomalous, logarithmic law |X n | — ao(|E|) ln(n). On comparing 
our results with earlier predictions by Brummelhuis and Hilhorst (J. Stat. Phys. 53, 249 
(1988)) for the TP diffusivity D n in the unbiased case, we infer that the Einstein relation 
fi n = /3D n between the TP diffusivity and the mobility p n = lmi| E |^ (|X n |/|Ejro) holds in 
the leading in n order, despite the fact that both D„ and fi n are not constant but vanish 
as n — > co. We also generalize our approach to the situation with very small but finite 
vacancy concentration p, in which case we find a ballistic-type law |X n ] = irao(\E\)pn. We 
demonstrate that here, again, both D n and fj. n , calculated in the linear in p approximation, 
do obey the Einstein relation. 



1 Introduction. 



Consider a square lattice of which each site except one is filled with a hard-core particle. The empty 
site is referred to as a "vacancy". The particles move randomly on the lattice, their random walks 
being constrained by the condition that each site can be at most singly occupied. More specifically, 
at each moment of time n = 1, 2, 3, . . . one particle selected with probability 1/4 among the four 
particles surrounding the vacancy will exchange its position with the vacancy. Suppose next that 
one selects one of the particles, "tags" it and follows its trajectory X n . Evidently, dynamics of 
the tagged - the tracer particle (TP) will be quite complicated, in contrast to the standard, by 
definition, lattice random walk executed by the vacancy: The TP can move only when encountered 
by the vacancy and its successive moves will be correlated, since the vacancy will always have a 
greater probability to return for its next encounter from the direction it has left than from a 
perpendicular or opposite direction. On the other hand, it is clear that on a two-dimensional 
lattice the TP will make infinitely long excursions as n — > oo even in the presence of a single 
vacancy, since its random walk is recursive in 2D and the vacancy is certain to encounter the 
tracer particle many times. A natural question is, of course, what are the statistical properties of 
the TP random walk, its mean-square displacement X 2 t from its initial position at time moment 
n, and the probability Pi' r ^(X) that at time n the TP appears at position X = {x\,X2)l 

The just described model, which represents, in fact, one of the simplest cases of the so-called 
"slaved diffusion processes", has been studied over the years in various guises, ranging from the 
"constrained dynamics" model of Palmer |IJ, vacancy- mediated bulk diffusion in metals and crys- 
tals (see, e.g. [§-0), frictional properties of dynamical percolative environments [^,^| or dynamics 
of impure atoms in close-packed surfaces of metal crystals, such as, e.g., a copper ||Io|~|l2]|. Brum- 
melhuis and Hilhorst jl3| were first to present an exact solution of this model in the lattice 
formulation. It has been shown that in the presence of a single vacancy the TP trajectories are 
remarkably confined; the mean-square displacement shows an unbounded growth, but it does grow 
only logarithmically with time, 

X ™ 1 TV asri ^°°> (!) 

which implies that the TP diffusivity D n , defined as 

_ XJ_ ln(re) 
71 ~ An ~ 47r(7r-l)n' [ ) 

is not constant but rather vanishes as time n progresses. 

Moreover, it has been found |l3| that at sufficiently large times Pn^ (X) converges to a limiting 



form as a function of the scaling variable n = |X|/ ^/ln(n). Still striking, this limiting distribution 
is not a Gaussian but a modified Bessel function Kq(j]), which signifies that the successive steps 
of the TP, although separated by long time intervals, are effectively correlated. These results have 
been subsequently reproduced by means of different analytical techniques in Refs. Jlif and ]l5|, ft6|. 
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Brummelhuis and Hilhorst have also generalized their analytical approach to the case of a very 
small but finite vacancy concentration p p7[ , in which case a conventional diffusive- type behavior 

p ft 

X ™ = 7 TT> P^ 1 ' n^oo, (3) 

(tt - 1J 

has been recovered. Note that Eq.(||) coincides with the earlier result of Nakazato and Kitahara 
H in the limit p <C 1 , and is well confirmed by numerical simulations P,|l8||. 

This paper is devoted to the following, rather fundamental to our point, problem: Suppose that 
we charge the tracer particle, (while the rest are kept neutral), and switch on an electric field E. 
In such a situation, the TP will have asymmetric hopping probabilities and in its exchanges with 
the vacancy, depending on the TP and vacancy relative orientation, the TP will have a preferency 
(or, on contrary, a reduction of the rate) for exchanging its position with the vacancy compared to 
other three neighboring particles. One might expect that in this case the TP mean displacement 
X„ will not be exactly equal to zero and might define the TP mobility as 

p n = hm H. (4) 
|E|— >o |E|rc 

Now, the question is whether the mobility p n , calculated from the TP mean displacement in 
the presence of an external electric field, and the diffusivity D n , Eq.(^), deduced from the TP 

mean-square displacement in the absence of the field, obey the generalized Einstein relation of the 
form 

H n =(3D n , (5) 

where [3 denotes the reciprocal temperature? 

Note that this question has been already addressed within the context of the TP diffusion in 
one-dimensional hard-core lattice gases with arbitrary finite vacancy concentration f?],|l9|, pp|. 
It has been found that Eq.(j|) holds not only for the TP diffusion in a ID hard-core gas on a 
finite lattice ||, but also for infinite ID lattices with non-conserved and conserved particles 
number ||, [19], M. Remarkably, in the latter case Eq.(||) holds for n sufficiently large despite the 
fact that both the TP mobility and the diffusivity are not constant as n — > oo but all vanish in 
proportion to 1/y/n [l9],^(J. On the other hand, it is well known that the Einstein relation 
is violated in some physical situations; for instance, it is not fulfilled for Sinai diffusion ]29[ ] or 
diffusion on percolation clusters, due to effects of strong temporal trapping in the dangling ends 
(see also Refs.|j27| and |2g] for some other examples). Hence, in principle, it is not a priori clear 
whether Eq. (||) should be valid for the model under study; here, the TP walk proceeds only due to 
encounters with a single vacancy, its mean-square displacement grows only logarithmically with 
time and the diffusivity follows much faster decay law in Eq.(||), compared to the D n ~ n -1 / 2 law 
obtained for the one-dimensional systems with finite vacancy concentrations. 

The paper is structured as follows: In section 2 we present more precise formulation of the 
problem and introduce basic notations. In Section 3 we discuss our general approach to computa- 
tion of the probability Pn r ' (X) of finding the TP at position X at time moment n and to evaluate 
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Pn (X) in the general form as a function of some return probabilities describing the random walk 
executed by the vacancy. The Section 4 is devoted to calculation of these return probabilities in 
the general case, as well as to the derivation of explicit expressions determining their asymptotical 
behavior. In Section 5, we present explicit asymptotical results for both the probability distribu- 



appropriate scaling variables, converges to a rather unusual limiting distribution. We also demon- 
strate here that the TP mobility, which is obtained in the present work in the leading in n order, 
and the TP diffusivity in the unbiased case, calculated earlier by Brummelhuis and Hilhorst [ fl3|| , 
do obey the Einstein relation. Further on, in Section 6 we extend our approach to the situation 
with very small but finite vacancy concentration and determine, in the leading in n order, the TP 
mobility. We show that also in this case the TP mobility and the TP diffusivity in the unbiased 
case do obey the Einstein relation, in the linear in p approximation and in the leading in n order. 
Finally, in Section 7, we conclude with a brief summary and discussion of our results. 



Consider a two-dimensional, infinite in both x\ and xi directions, square lattice every site of which 
except one (a vacancy) is filled by identical hard-core particles (see Fig.l). All particles except one 
are electrically neutral. The charged particle, which is initially at the origin, will be referred to in 
what follows as the tracer particle - the TP. Its position at the lattice at time n will be denoted 
by X„. Electric field E of strength E = |E| is oriented in the positive x\ direction. For simplicity, 
the charge of the TP is set equal to unity. 

Next, we suppose that at each tick of the clock, n = 1,2,3,..., each particle selects at random 
a jump direction and attempts to hop onto the target site. Evidently, the jump event can be only 
successful for four particles adjacent to the vacancy. 

The form of the jump direction probabilities depends on whether the particle is charged or not. 
For uncharged particles all hopping directions are equally probable and hence, all jump direction 
probabilities are equal to 1/4. On the other hand, the charged particle - the TP, "prefers" to jump 
in the direction of the applied electric field; the normalized jump direction probabilities of the TP 
are given, in a usual fashion, by 



where Z is the normalization constant, e v is the unit vector denoting the jump direction, v G 
{±1,±2}, and (E • e^) stands for the scalar product. We adopt the notations e±x — (±1,0) 
and e± 2 — (0,±1), which means that e 1 (e_x) is the unit vector in the positive (negative) X\- 
dircction, while e% (e_2) is the unit vector in the positive (negative) 22-direction. Consequently, 
the normalization constant Z is 



tion and the TP mean displacement. We show that as n — ► 00 



, (X) written in terms of two 



2 The model. 



p, y = Z 1 exp -(E ■ e v ) , 



(6) 




(7) 
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where the sum with the subscript /i denotes summation over all possible orientations of the vector 
e M ; that is, fi = {±1,±2}. Note that the jump direction probabilities defined by Eqs.Q and (|^) 
do preserve the detailed balance condition. 

E 

► 




Figure 1: Two-dimensional, infinite in both directions, square lattice in which all sites except one are 
filled with identical hard-core particles (grey spheres). The black sphere denotes a single tracer particle, 
which is subject to external field E, oriented in the positive xi direction, and thus has asymmetric hopping 
propabilities. 



Next, it is expedient to reformulate the dynamics between two consequtive jumps of the TP in 
terms of the random walk executed by the vacancy and its jump direction probabilities q v . From 
the viewpoint of the vacancy, the jump direction probabilities depend on whether the TP is one 
of four surrounding particles or not. Evidently, in case when all four surrounding particles are 
electrically neutral, we still have that on the next time step the vacancy will change its position 
with one of four neighboring particles selected at random with equal probabilities. Hence, in case 
when the TP is not adjacent to the vacancy, all four jump directions for the vacancy are equally 
probable, i.e. q v = 1/4. On the other hand, the situation is a bit more complex when one of these 
four particles is the TP, which has asymmetric jump direction probabilities, Eq.(^|). A natural 
choice of the normalize^ jump direction probabilities of the vacancy in this case is as follows: 
Suppose that at time moment n the tracer particle is at position X„ and the vacancy occupies 

1 Note, that normalization here insures that the vacancy performs one jump each time step. Otherwise, we will 
introduce artificial "temporal trapping" probability, which would definitely lead to the violation of Eq. (ph. 
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an adjacent site X„ + e„. Then, an exchange of the positions between the TP and the vacancy, 
which implies that the TP is moved one step in the e„-direction, takes place with the probability 

q_ v = Z* p v (8) 

while the probability of the exchange of positions with any of other three adjacent particles is 
given by 

%±-v = \z* (9) 

The normalization constant Z* in this case is, evidently, 

Z*=3/4 + p„, (10) 

where p u has been defined previously in Eq.(||). 

Consequently, apart of four sites in the immediate vicintiy of the tracer particle, the vacancy 
performs a standard, symmetric random walk. In the vicinity of the TP, the vacancy jump 
direction probabilities are perturbed by the TP asymmetric hopping rules. Hence, the random 
walk executed by the vacancy can be thought off as a particular case of the so-called "random 
walk with defective sites" (see Ref.|| for more details), or as a realization of the "random walk 
with a hop-over site" jl6| . 



3 Probability distribution function 

A standard approach to define the properties of the TP random walk would be to start with a 
master equation determining the evolution of the whole configuration of particles. In doing so, 
similarly to the analysis of the tracer diffusion on 2D lattices in the presence of a finite vacancy 
concentration (see, e.g. Ref.^6|), one obtains the evolution of the joint distribution P n (X, Y) of 
the TP position X and of the vacancy position Y at time moment n. The property of interest, i.e. 
the reduced distribution function of the TP alone will then be found from P„(X, Y) by performing 
lattice summation over all possible values of the variable Y. 

Here we pursue, however, a different approach, which has been first put forward in the original 
work of Brummelhuis and Hilhorst pUf ; that is, we construct the distribution function of the TP 
position at time n directly in terms of the return probabilities of the random walk performed by 
the vacancy. The only complication, compared to the unbiased case considered by Brummelhuis 
and Hilhorst p3| , is that in our case ten different return probabilities would be involved, instead 
of three different ones appearing in the unbiased case. Hence, the analysis will be slightly more 
involved. 

We begin by introducing some basic notations. Let 

• Pn^ (X) be the probability that the TP, which starts its random walk at the origin, appears 
at the site X at time moment n, given that the vacancy is initially at site Yq. 
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• F*(0 | Yo) be the probability that the vacancy, which starts its random walk at the site Yo, 
arrives at the origin for the first time at the time step n. 

• F*(0 | e„ | Y ) be the conditional probability that the vacancy, which starts its random 
walk at the site Yo, appears at the origin for the first time at the time step n, being at time 
moment n — 1 at the site e„ . 

Further on, for any time-dependent quantity L n we define the generating function of the form: 

+00 



L(0 = J2 L nC (11) 



71=0 

and for any space-dependent quantity Y(X) the discrete Fourier transform 

y(k)=^exp(i(k-X))F(X), (12) 
x 

where the sum runs over all lattice sites. 

Now, following Brummelhuis and Hilhorst Jl3| , we write down directly the equation obeyed by 
the reduced probability distribution P„ (X) (cf Ref.[jl5| for a study of the joint probability of 
the TP position and of the vacancy position in the unbiased case): 



P(*-)(X) = 5 x ,o 1-X>/(0|Y ) 

+ X X ■■■ X X <5 mi+ ... +mp+1 ,„X---X (5e -i+-+ e - P ' x x 

p—l mi — 1 m p — 1 m p +i— v\ v v 

i- X ^(°i- e - P ) 

3=0 

x K»(p\ev p | -e„ p _J...i£ a (0|e I , a | -cOF'^Ole^ | Y ). (13) 
Next, using the definition of the generating functions and of the discrete Fourier transforms, 



Eqs.([l l|) and (12), we obtain the following matricial representation of the generating function of 
the TP probability distribution: 



(14) 



P^)(k;0 = ^ ^1 +P- 1 (k;e)E t/ M(k;^*(0 | e p | Y ;o) ■ 
In Eq.(|T4|) the function 2?(k; £) stands for the determinant of the following 4x4 matrix, 

2>(k;0=det(I-T(k;0), (15) 
where the matrix T(k; £) has the elements ( T(k; £) J defined by 

(T(k;0) v = exp(i(k-e„))i4„,_ p (£). (16) 



G 



T(k;0 



(17) 



Explicitly, the matrix T(k; f) is given by 

I e^ViK) e*4y(0 e'^yag) e w Mi i2 (0 \ 

3-* fc M_ lf _i(^) e- <fcl A_i,i(0 e-**M_i,_ 2 (0 e- <fcl A_i i2 (0 

e^iVitf) e*"A a ,i(0 e* fc M 2 ,_ 2 (£) e lfc M 2 , 2 (£) 

Ve-^A-a,.!^ e~^A_ 2A (C) e- 4fc M_ 2 _ 2 (£) e~^A^ 2 (Oj 

where the coefficients A v ^{^)^ u,fx = ±1,±2, stand for 

+00 

4,,„(0 = F*(0 I e„ I e M ;0 = £)F*(0 | e„ | e M )f , (18) 

n=0 

i.e. are the generating functions of the conditional probabilities for the first time visit of the origin 
by the vacancy, conditioned by constraint of the passage through a specified site on the previous 
step. Note that, by symmetry, 

A2AO =A-2,A0, 

il„,2(0= 4,-2(0 (19) 

for v — ±1 and 

Aa )2 (0=A_ 2 ,_ 2 (0, 

4.-2(0=^-2,2(0- (20) 

As a result of such a symmetry, we have to consider just ten independent functions A^ ;1/ (0 (note 
that in the unbiased case one has to deal with only three such functions |lq] ). Explicit expression 
of the determinant in Eq. (|l5|) in terms of these generating function is presented in the Appendix. 
Lastly, the matrix C/ M (k; £) in Eq.([l4|) is given by 

C/^k;0 = 2?(k;0E( 1 - ^ k ^))(I - T(k;0)^ e 1 ^. (21) 

The property of interest - the TP probability distribution function, will be then obtained by 
inverting p(* r )(k; with respect to the wave- vector k and to the variable £ : 



(22) 



where the contour of integration C encircles the origin counterclockwise. 

Finally, we remark that as far as we are interested in the leading large-n behavior of the prob- 
ability distribution Pn^ (X) only, we may constrain ourselves here to the study of the asymptotic 
behavior of the generating function P(* r )(k; £) in the vicinity of its singular point nearest to £ = 0. 
We notice that similarly to the unbiased case, this point is £ = 1 when k = 0. As a matter of 
fact, such a behavior stems from the a priori non-evident fact that the vacancy, starting from 
a given neighbouring site to the origin, is certain to eventually reach the origin. This will be 
demonstrated explicitly in section 4 (cf. Eq.jrij)); as a matter of fact, one can see from Eq.(|47|) 
and the explicit representation of 2?(0;0 presented in the Appendix that 2?(0;£ = 1) = 0. In 
consequence, expansion in powers of a small deviation (1 — £) has to be accompanied by a small-k 
expansion, exactly as it has been performed in Ref. Il3|. 
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if s' 


i {0,e ± i,e ±2 }, 


5 s ,o - 1/4 


if s' 


= 0, 


5q v 


if s' 


= e„ and s = 0, 


-<5<W 3 


if s' 


= e„ and s' 7^ 0, 



4 The return probabilities F*(0 | e M | e„). 

As we have already remarked, the vacancy random walk between two successive visits of the lattice 
site occupied by the TP can be viewed as a standard, two-dimensional, symmetric random walk 
with some boundary conditions imposed on the four sites adjacent to the site occupied by the TP. 
In order to compute the return probabilities F*(0 | e M | e„) for such a random walk, we add, in 
a usual fashion J5|^], an additional constraint that the site at the lattice origin is in absorbing 
state. Then, the vacancy random walk can be formally represented as a lattice random walk with 
site-dependent probabilities of the form p + (s|s') = 1/4 + g(s|s'), where s is the site occupied by 
the vacancy at the time moment n, s' denotes the target, nearest-neighboring to s site, 



q( S \s')={ ■ (23) 



where 5q v is defined, according to Eqs.(||),(||) and (|Io|), by 

^-^m~\ « 24 > 

Further on, we define P„(s | So) as the probability distribution associated with such a random 
walk starting at site So at step n = 0. 

Now, let the symbols £, A and B define the following three events: 

• the event £: the vacancy, which has started its random walk at the site e„, visits the origin 
for the first time at the n-th step exactly, being at the site e M at the previous step n — 1; 

• the event A: the vacancy, which started its random walk at the site e„, is at the site e M at 
the time moment n — 1 and the origin has not been visited during the n — 1 first steps of 
its walk; 

• the event B: the vacancy jumps from the neighboring to the origin site e M to the site at 
the n-th step exactly. 

Evidently, by definition, the desired first visit probability F^(0 | e M | e^) is just the probability 
of the £ event 

F*(0|e M |e v )=Prob(£). (25) 
To calculate Prob(£ ) we note first that the probabilities of such three events obey: 

Prob(£) = Prob(^ nB) = Prob(^) Prob(S). (26) 
On the other hand, we have that 

Prob(^l) = P+_ x (e„ I e„), (27) 
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and 

Pr ° b(S) = (28) 

Hence, in virtue of Eqs.(p5[),(p^),(p7|) and (28), the return probability F*(0 | e M | e„) is given 
explicitly by 

F*(0 I e M I e„;0 = £ (jjff^) ^ + ( e M I e„;£)- (29) 

Therefore, calculation of the return probabilities F*(0 | e M | e„) amounts to the evaluation of the 
probability distribution P^~(s | So) of the vacancy random walk in the presence of an absorbing 
site placed at the lattice origin. Such a probability distribution will be determined in the next 
subsection. 

4.1 The generating function of the probability distribution F + (s | s ). 

Making use of the generating function technique adapted to random walks on lattices with defective 
sites || and we obtain 

2 

P+( Si I sj; = |sj;0+ J2 A ^ I "I? ^) F+ ( Sl I S J'^)' ( 3 °) 

l=-2 

where 

f ei , forie{±l,±2}, 
Si = < (31) 

0, for i = 0, 

and 

A( Si | Sl ;C)=^ P ( Si l s ';09(s'|si), (32) 

s' 

P(si I sj;£) being the generating function of the unperturbed associated random walk (that is, 
symmetric random walk with no defective sites). 

Further on, Eq. ( p0| ) can be recast into the following matricial form: 

P+ = (1-A)- 1 P, (33) 

in which equation P, P + , A stand for the 5x5 matrices with the elements defined by 

P <ii = P(s i |s j ;0, P + *j = ^ + (si I sj;0, A <ii = A( 8l |8 J ;0 > (34) 
where i,j — 0, +1, —1, +2, —2. Using next an evident relation ||: 

p(s k 1 si; - s k j + p ( Sk 1 si + e -; 0. (35) 

V 

and the symmetry properties of a standard random walk, one can readily show that: 



9 



for si ^ s and Sk ^ s , 



A(s k | s i; = 3 Hi ( P(0 | 0; - 1 - | s i; + ) , 



for si ^ Sq and Sk = So, 



1 



for si = s , 



A(s | bi; = ( P(0 | 0;0 — |a ( P (° I ~ 1 



A(s k | s ; - 4,o - (1 - 0P(sk I 0; 0. 



(36) 



(37) 



(38) 



Consequently, the matrices A and P in Eq.(|33|) are given by 

( a 5qif Sq-xf Sq 2 f Sq 2 f \ 

b Sq—ie Sq2C 5q 2 c 

b 5q%e Sq 2 c Sq 2 c 

b 8q±c Sq-ic 5q 2 e 

\b Sq\C 8q-\C Sq 2 e 0, 



where 



a ee 1-(1-0G(0, b=±j£(l-G(Q), e ee 1(2^(0-1), 
c = |(-l + ^ + 2G(0(l-^)-fl(0), 



and 



with 



' G(0 


(G(o-i)A (G(o-i)/e 


- 1)A 


- 1)/A 


(o(o - i)/e 








- 


G(0 - 2g(0 G(0 




r(0 


(G(0 - 




G(0 




\(G(o - i)/e 








G(0 = 


P(0|0;0, 9(0 = ^(^(ei 


- ei ;£)-P(0 


|0;0). 



r(0 



(39) 



(40) 



(41) 



(42) 



Note that Eqs.(^9|) and ([i"l| ) now define the P + matrix explicitly, and hence, define the generating 
function of the probability distribution P + (s | sq). 



4.2 Asymptotic behavior of the generating functions of the return prob- 
abilities in the vicinity of £ = 1. 

As we have already remarked, here we constrain our consideration to the analysis of the leading 
in n behavior; this amounts to consideration of the leading in the limit £ — > l~ behavior of the 



10 



corresponding generating functions. Expanding G(£) and g(£) in the vicinity of the singular point 
£ = 1, (cf Refs.g and |D|||,|j|), we have 

G(0 = iln-^--i(l-£)ln(l-£) + O(l-£) ! £-1-, (43) 

7T 1 — £ 27T 

and 

5 (£) = f2 - -) + -(1 - ln(l -0+O ((1-0), £ - I". (44) 

Consequently, we find by solving the matricial equation (|3^) , that the generating functions of the 
return probabilities obey 

• 4 -«»-^ir-w( 1 "< 1 -«)" +0 < 1 -«- < 45 > 

where u = exp(fiE/2), A$,(u) and A%l(u) are some rational fractions (all listed explicitly in the 
Appendix), while 

S(u) = I (7T - 2) M 6 + (2 7T 2 - 6 7T + 12) U 5 + (8 7T 2 - 25 7T + 34) u 4 - 

- (4tt 2 - 60 tt + 88) u 3 + (8tt 2 - 25tt + 34)u 2 + (2tt 2 -6tt+ 12^u + 7T- 2}. (46) 

It follows from Eqs.([45|) and explicit expressions for A„ ift (£) presented in the Appendix, that, in 
particular, the generating functions of the return probabilities fulfil: 

Ai,_i(l-)+i4_i,_i(l-) + 2A 2l _i(l-) = 1 
A 1)1 (l-)+A_ 1)1 (l-) + 2A 2 , 1 (l-) = 1 

A li 2(l-)+A-i,2(l-)+A_ 2i2 (l-)+i42,2(l-) = 1, (47) 

which relations imply that the vacancy, starting its random walk from a given, neighbouring to 
the origin site, is certain to return eventually to the origin. 

5 The TP mean displacement and the probability distribu- 
tion. 

In this section, we proceed as follows: Taking advantage of the asymptotical expansion obtained 
in the previous section, we first determine the small (1 — £) behavior of the generating function 
p(* r )(k;£), accompanied by the small-k expansion. Next, we evaluate the generating function of 
the TP mean displacement, by differentiating the obtained asymptotical expression for p( tr ) (k; £) 
with respect to the components of the wave-vector, and analyse its large-n behavior. Lastly, we 
invert the asymptotical expansion of the generating function p( tr )(k; £) and obtain the correspond- 
ing probability distribution P„ r ^(X) in a certain scaling limit. 
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5.1 Asymptotic expansion of the generating function P( tr )(k; £) . 

Using the explicit representation of the determinant 2?(k; £) in Eq.(|T5|) in terms of the generating 
functions of the return probabilities A;, iAt (£), presented in the Appendix, as well as the asymptotical 
expansions in Eq. ( p5| ) , we find that in the vicinity of £ = 1 and for small values of the wave- vector 
k, X>(k;£) is given by 

2?(k; = i^i(u)fci + T 2 {u)k\ + F 3 (u)k% - Fa{u) In" 1 (1 - f) + . • • , (48) 

where we have used the shortenings 

^ {U) = ( U 2 + 2(7r-l)u + l)5(u) ' (49) 

(tt - 2)(1 + M ) 4 (» 2 + 1)(?/ 2 + 2(2^ - 3)u + 1) 

^ « = of 2 , or ^ i i\ar \ ' ( 50 ) 

2{u z + 2{n — l)u + 1)6 (u) 

u(tt - 2)(1 + u)\(2tt - j)u 2 + 2 u + 2 tt - 3) 

^3(«) = 7-0 , „, 7T- , 7T^7-T > ( 51 ) 



and 



(ir-2)(u- 1)(1 + u) 5 (u 2 - 


f 2(2tt- 


3)« + l) 


(w 2 + 2(tt-1)w- 


f- l)S(u) 




(tt-2)(1 + u) 4 (u 2 + 1)(u 2 - 


f 2(2tt - 


3)u+l) 


2(u 2 + 2(tt-1)wH 


- l)S(u) 




u(tt- 2)(1 + m) 4 ((2tt- 3)u 


2 + 2u + 


2vr- 3) 


(u 2 + 2(?r - l)w + 


l)S(u) 




+ m) 4 ((2tt-3)m 2 + 2u + 2 


7T — 3)(l< 


2 + 2(2tt - 


(m 2 + 2(tt-1)m + 







Mu) = ' " ,, 2 ' i ' > '- '-, (52) 



and assumed, for simplicity, that the starting point Yo of the vacancy random walk is Yo = e_i. 
On the other hand, we find that 

J2u„(*0F*(0\e» I -e 1 ;0 = -^i(*-^-^(#2+-- ■ (53) 

Consequently, in the small-k limit and £ — ► 1~, the generating function P(* r )(k;£) obeys 
P^(k; = ^{l - (-ia h + ^ ai kf + \a 2 k^j In (1 - 



(54) 



where the coefficients 

a (E) = tt- 1 smh(/3E/2)((2n - 3) cosh(/3£/2) + 

ai (E) ee tt- 1 cosh(/3S/2)((27r - 3) cosh(/3£/2) + l)" 1 , (55) 
a 2 {E) ee 7r- 1 (cosh(/3£;/2) + 2tt - 3)"\ 
are all functions of the field strength E and of the temperature only. 

5.2 The TP mean displacement for arbitrary field strength E. 

As a matter of fact, the leading large- n asymptotical behavior of the TP mean displacement can 
be obtained directly from Eq.(|54|), since the generating function of the TP mean displacement, i.e. 

n=0 
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obeys (see, e.g. Ref.Q): 

^)=-^(0;e)e 1 + ^(0;e)e 2 j. (57) 

Consequently, differentiating the expression on the right-hand-side side of Eq. ( |54| ) with respect to 
the components of the wave-vector k, we find that the asymptotical behavior of the generating 
function of the TP mean displacement in the vicinity of £ = l - follows 

X(0~(^lnAK (58) 



Further on, using the discrete Tauberian theorem (cf. Ref. Q) and Eq. (|55|) , we find the following 
general force-velocity relation for the system under study 



X 



/ , \ fl sinh(/3£72) , \ 

» ~ (a (E) Inn) ei = ( ~ (2?r _ 3) CQsh ( ^ /2) + ! H as n ^ oo, (59) 



which shows that the TP mean displacement grows logarithmically with n. In consequence, one 
may claim that the typical displacement along the x± direction scales as ln(n) as n — > oo. On 
the other hand, typical displacement in the 22-direction is expected to grow only in proportion to 



y/hi{n) , as in the unbiased case |13|. These claims will be confirmed in what follows by the form 
of the scaling variables involved in the limiting distribution. 

Consider next behavior of the coefficient <Xq(E) in the limit E — ► 0. Here, we find from Eq.(|55|) 
that 

<xo(E) = j f E U +0(E 3 ), (60) 

and hence, the mobility defined in Eq.(|J), follows 

P ln ( n ) fats 

fJ>n ~ ~. — 7T , as n ► oo (61) 

47r(7r — 1) n 

Comparing next the result in Eq.(|6l|) with that for the diffusivity D n , Eq.(||), derived by Brum- 
melhuis and Hilhorst in the unbiased case, we infer that the TP mobility and diffusivity do 
obey, at least in the leading in n order, the generalized Einstein relation of the form \x n = (3D n 
||. Note, that this can not be, of course, an a priori expected result, in view of an intricate na- 
ture of the random walks involved and anomalous, logarithmic confinement of the random walk 
trajectories. 

5.3 Probability distribution Pi ir) (X). 

We turn next to calculation of the asymptotic forms of the probability distribution Pn (X). 
Inverting p( tr ^(k;£) with respect to k, we notice first that in the limit £ — * 1~ the integrand is 
sharply peaked around k = 0, such that the bulk contribution to the integral comes from the 
values ki — and ki = 0; this implies that we can extend the limits of integration from ±7r to 
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±oo, which yields in the limit £ — ► 1 : 

f+oo r-\-oo 



(62) 



(63) 



P (tr) (X;£) ~ - / / dfcidfc 2 exp ( - ik\X\ - ik 2 x 2 ) x 

(1 -0(2tt) 2 J_ oc J_ ocl V / 

x {l- (-ia {E)k 1 + ^ai(E)kl + ^a 2 (E)k^ In (1 - 
Further on, using the integral equality 

{l - (-iaoOE)fci + ^ai(£)fc 2 + la 2 (E)k£j In (1 - o}^ = 

J dv exp (- «{l - f-<oo(S)*i + + ^2(^2) In (1 - £)}) 

we cast the integral in Eq.(|62]) into the form: 

P (tr) (X;£) - / dvexp(-u) / dfci / dfc 2 exp(-i/c 2 x 2 ) x 

x ex P (|(ai(£;)fc 2 + a 2 (£;)fc|)ln(l-e)-ifci(a: 1 +va (£;)ln(l-0)) (64) 

Note now that in order to evaluate explicitly the Gaussian integral in Eq.(p4[), we have to consider 
separately two cases: when (a) the external field in infinitely strong, E — oo (which implies 
a 2 = 0), such that the TP performs a totally directed walk, and (b) - when E is bounded, E < oo 
(and hence, a 2 > 0). 

5.3.1 Directed walk, E = oo. 

We start with the simplest case when the TP performs a totally directed walk under the influence 
of an infinitely strong field. In this case, the probability distribution is defined for non negative 
xi values only, and the equation ( |64| ) reduces to: 

6(x 2 )9(x 1 ) r+oo r+oo 



x 



exp {^ ai {E)kl In (1 - f) - iki (xt + va Q (E) In (1 - £)) ) (65) 



where 8{x\) denotes the Heaviside theta-function. Performing the integrals, we find that, in the 
limit £ — > 1~, the generating function p(* r )(X;£) obeys: 

Applying next the discrete Tauberian theorem we find eventually, 

e >(X,~ W 0^exp(-«4 (67, 

which means that in the totally directed case, in the large-n and large- a; i limit, the scaled variable 
rjoo = tt(2tt — 3)xi/ln(n) is asymptotically distributed according to 

■P(»?oo) = 0(Voa) exp(-Voc), (68) 
i.e. has an exponential scaling function. 
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5.3.2 Arbitrary bounded field E < oo. 

In this case the coefficient > and the probability distribution is defined also for negative 
values of x\ ; as well, P(' r )(X;£) is defined also for non-zero values of x-i- In this general case, we 
find, performing integrations over the components of the wave-vector, that p( tr '(X;£) attains, as 
£ — > 1 — , the following form: 

f +00 



x exp 



p( tr )(X; ~ -(2tt(1 - £) In (1 - V^P)«^)J / d«exp(- V )x 

The integral in the latter equation can be calculated exactly, which yields 

PW(X; ~ - (tt(1 - In (1 - Ov^M^))"* exp (^jf^i) #0 (tV)) < ( 70 ) 
where Kq is the modified Bessel (McDonald) function of zeroth order, and 



Finally, using the discrete Tauberian theorem jE], |TJ , we find from Eq. ( fzT| ) that in the large-n and 
large-X limits, the probability distribution P„ ^(X) obeys 

Note that in the unbiased case, i.e. when E = 0, the probability distribution P„ (X) defined by 



Eqs.(|7l|) and (|72| ) reduces to the form predicted earlier by Brummelhuis and Hilhorst [13 



5.3.3 Limiting probability distribution function. 



Now, we recollect that the scaling behavior expected is x\ ~ In (71) (for E > 0)and X2 ~ yni(n). 
In order to obtain from Eqs.(|7l|) and ( f72|) the limiting probability distribution, we introduce two 
scaling variables: 

r] 1 =x 1 /a Q {E)\n{n), 
m = X2/ v / 2ct2(E) In (n). 

Note that 771 becomes rjac in the special case E = 00. In terms of these scaling variables T]e(ti) in 
Eq.(^TJ) takes the form: 

ai(A) V ao(£/)ln(n)\ 771 / / 

Note now that for arbitrary fixed 771 and 772, the argument of the Bessel function 77^(71) written 
in terms of the scaling variables tends to infinity as n — > 00. Consequently, using the limiting 
behavior of the modified Bessel function 

1 /2 

Ko(2/)= © eM-y)(± + 0(l/y)), (75) 
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we find that the probability distribution P„ '(X), written in terms of the scaling variables, con- 
verges as n — > oo to the limiting form 

J (2TTa 2 (E)a 2 (E) m \n 3 (n)r 1 / 2 exp(~ m -r 1 2 2 / m ), for m > 0, 
for 771 < 0, 



P(, 1 ,rt=^ exp (_, 1 _^) (77) 



or, equivalently, that the scaling variables 771 and 772 have the following, rather unusual limiting 
joint distribution function: 

6>(?7l) / 77^ 

— = exp - 7/1 

/7T7yi v 771- 

We note that this distribution is properly normalized and yields, of course, the same result for the 
TP mean displacement as the approach based on differentiation of the asymptotical expansion of 
the generating function. We also remark that the reduced distributions P(r]i) — J drfyP{r)i,r}2) 
and P(rj2) — J dri\P(ri\, 772) take the form: 

P( m ) = exp(-2|77 2 |), (78) 
and hence, the reduced distribution P(rji) appears to be exactly the same as in the case E = 00. 

6 Finite vacancy concentration. 

In this last section, we generalize our analysis of the biased TP mean displacement to the case 
when vacancies are present at a very small, but finite concentration p. In our approach, we follow 
closely that of Brummelhuis and Hilhorst |l7j , who pointed out that for the unbiased case in the 
limit of low vacancy concentration, the many- vacancy problem can be interpreted in terms of the 
one-vacancy solution, which entails meaningful results to the leading order in the concentration 
of vacancies for p <C 1. We thus just extend here their consideration over the biased case. 

Following Ref. JlTj , we begin by considering a finite lattice of size LxL, containing M vacancies. 
The mean concentration of the vacancies is thus p = M/L 2 <C 1. We suppose that the charged TP 
is initially at the origin and initial positions of the vacancies are Yg 1 '' , Yq 2 ' , . . . , Yq A/ ' 1 , which all 
are different from each other and from 0. All other sites are filled with neutral hard-core particles. 
The field E is again supposed to be oriented in the positive xi-direction. 

Similar to the single vacancy case, we stipulate that at each time step, all vacancies exchange 
their positions with either of neighboring particles, such that each vacancy makes a step each time 
step. Exchanges with the charged TP are governed by the same rules as described in Section 2. 
Note that, of course, when many vacancies are present, it may appear that two or more vacancies 
occupy adjacent sites or have common neighboring particles, in which case their random walks will 



interfere. However, as noticed in |17|, these cases contribute only to 0(p 2 ) and thus go beyond 



our approximation; hence, we discard such possibilities here. 
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Now, let 7 , i* r ' ) (X|Yg 1 ' ) , Yq 2 '', . . . , Yq M ^) denote the probability of finding at time moment n 
the TP at position X as a result of its interaction with all M vacancies collectively. Further on, 
let Pn (X^|Yq ) denote the probability of finding the TP at site X^ at time moment n due to 
interactions with a vacancy initially at Yq in a system with a single vacancy. Then, assuming 
that the vacancies contribute independently to the TP displacement, one has, following Rcf . |jl7f : 

M 

^(X|Y« y« , . . . , Y^) «£■•• E <W + ... + v<"> II ^ r) (X«|Y«) (79) 

1 1 

Upon averaging over all initial vacancy configurations, and denoting this average by the angular 
brackets, we have 

M 

(^>cx|y« yw ... ,y< m) )) « E ••• E ^ + ... +Y r<n^ tr) (x«iYW ) ) ! (80) 

1 1 

which, in the limit of very small vacancy concentration, simplifies to |l7j ]: 

M 

(p<^(X|Y« Y?», . . . , YD) - E • • • E <W> + ... +Y -> II (^ tr) (X«|Y«)) (81) 

V (l) V (M) j = l 

1 1 

Now, defining the Fourier transformed distributions 

T>W(k,M,L) = Eexp(z(k • X))(^ ) (X|Y«, Y< 2) , . . . , Y<, M) )), (82) 
x 

and 



Pf-)(k) =Eex P (z(k.X))P r f)(X|Y«), (83) 



and performing the corresponding summations, we find that 

£<M(k,M,L)« (pi tr) (k)) M (84) 
Turning next to the limit L, M — * oo (while the ratio M/L 2 = p is kept fixed), we obtain 

(k, p) = Jim P r f > (k, M, L) « exp ( - p n„ (k)) , (85) 

where 

n 

n n (k) = EE A ^( k i e ") E (86) 

A„(k| ei/ ) = l-P f C* r )(k|-e v )exp(t(k.e tf )) 1 (87) 

F*(0|e„|Y) are conditional return probabilities, defined in Section 3, and Pi* r ' ) (k|— e^) is the 
Fourier transformed single- vacancy probability distribution P„ (X|— e„). The latter can be read- 
ily obtained by applying the discrete Fourier transformation to the relation 

P^(X|Y) = 5 x ,o fl-E*7(°l Y ) 

n 

3=0 v 
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and chosing Y = — e u . 

Further on, using the results of the previous section, we find that in the limit £ — > 1 _ , k — > 0, 
the generating function of O n (k) is given by 



with 



n(k;0 = A(k|e v ;0 £ ^*(0|e„|Y; £), 

v Y^O 



A(k| ei/ ;0 = T ^(l-e X p(i(k-e I/ )){l-ln(l-0(-ia (E)A;i 
+ i ai (£;)fc? + ia 2 (S)fe|)} _1 ) +... . 



(89) 



(90) 



We turn next to calculation of ^2y-t F* (0\e„\Y; in the limit £ — ► 1 , k — * 0, which can be 
done rather straightforwardly by taking advantage of the results of Section 3. We have then 

Vv 



£F*(0| ei/ |Y;£)=£ 

Y#0 



3/4 + p 



3/4 + Pv 
-) Bta-AJ-^BCYjO, 



Y^O 



(91) 



Y/0 



where is the v — th basis vector, B* denotes the transposition of 2?„, and S(Y; £) is the vector, 
whose elements are (P(si|Y; i = 0,1,-1,2,-2. Explicitly, 
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Further on, using an evident symmetry relation 

P( fli |Y i ;0=i > (Y i |«,;0, 



(92) 



as well as the relation in Eq.(p5|), we obtain 
1 „,.A „ ( 1 



Y^O 



G(0 Bo 



-(G(0-1) (B1 + B-1+B2 + B-2) 



(93) 



(94) 



Then, combining Eqs.(J9l|), (|94j) and (|39j), (J40]), (f43|), and performing some straightforward but 
cumbersome calculations, we find that in the limit £ — > l - and k — * 0, the sum X^y^o F* (0|e„ | Y; £) 
is given by 

7T 



^F*(0|e v |Y;£) = -- 



Y^O 



(i-OMi-0 



(95) 



which is, remarkably, independent of u and in the leading in £ order. Consequently, in the limit 
£ — > 1~ and k — * 0, the generating function f2(k;£) obeys: 

7T -ia (E)ki + \on{E)kl + ±a 2 {E)k% 



n(k;f) 



(1 - £) 2 1 - In (1 - {-iao{E)k x + ±ai(P)fc 2 + |a 2 



(96) 
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Next, using the discrete Tauberian theorem, we obtain from the latter equation that in the limit 
n — > oo and k — > 0, 



( - iao{E)k! + ^ai(E)kf + \a 2 {E)k 



n «( k ) %7r 1 + m ( n )(-ia (E)k 1 + la 1 (E)k 2 1 + ^a 2 (E)kl) (9?) 



Finally, inverting Eq.(|8^) with respect to the wave- vector, 

^ tr) (X,p)«^y dkt J dfc 2 exp(-i(k-X)-pn n (k)), 



(98) 



and taking advantage of Eq.(p7|), we find that the leading, large- n behavior of the TP mean 
displacement is given by 



X,i 



(n ao (E) P n) ei = _ + e,, (99) 



i.e. grows linearly with time. This signifies that the TP mobility attains a constant value at 
sufficiently large times n, 

Mn = Um 1^1 = ^ (100) 
P |BHo |E|n 4(tt-1) v ' 

Lastly, on comparing Eq. ( |100| ) and the result of Brummelhuis and Hilhorst |l7j] for the TP diffu- 
sivity in absence of the field, Eq.(pl), we notice that again the Einstein relation is fulfilled! 



7 Conclusion 

In conclusion, we have studied the dynamics of a charged tracer particle diffusing on a two- 
dimensional lattice, all sites of which except one (a vacancy) are filled with identical neutral, 
hard-core particles. The system evolves in discrete time n, n = 0, 1, 2, . . . , by particles exchanging 
their positions with the vacancy, subject to the condition that each site can be at most singly 
occupied. The charged TP experiences a bias due to external field E, which favors its jumps 
in the preferential direction. We determine exactly, for arbitrary strength of the field E — |E|, 
the leading large-n behavior of the TP mean displacement X„, which is not zero here due to 
external bias, and the limiting probability distribution of the TP position. We have shown that 
the TP trajectories are anomalously confined and its mean displacement grows with time only 
logarithmically, X„ = (ao(E) ln(n)) ei as n — * oo. On comparing our results with the earlier 
analysis of the TP diffusivity D n in the unbiased case by Brummelhuis and Hilhorst |D| , we have 
demonstrated that, remarkably, the Einstein relation \i n = (3D n between the diffusivity and the 
mobility \x n of the TP holds in the leading in n order, despite the fact that both D n and fj, n tend 
to zero asn^ oo. Note, however, that validity of the Einstein relation for the system under study 
relies heavily on the proper normalization of the vacancy transition probabilities (see, Eqs.(||),(|^) 
and (|ic|)). In absence of such a normalization, artificial "temporal trapping" effects may emerge, 
which will result ultimately in the violation of the Einstein relation for the system under study (see 
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also Refs.[g7| and |£8[ for physical situations in which such type of effects is observed). Further on, 
we have also generalized our approach to the situation with small but finite vacancy concentration 
p, in which case we have found a ballistic-type law of the form X„ = (tt ao(E) p n) e\. We have 
shown that here, again, both D n and p, n calculated in the linear in p approximation do obey the 
Einstein relation. 
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8 Appendix. 

In this Appendix, we list some explicit expressions skipped in the body of the manuscript. First 
of all, explicit form of the determinant X>(k;£) in Eq.(j^) in terms of the generating functions of 
the return probabilities A v ^ = -4„, M (£) reads 



2?(k;0 = 1 - A 2 2 +2A_i,_iA 2) iA 1) 2A_a,2 + 2A 1 ^ 1 A 2il A_ lt2 A 2t2 - 

- 2 A 2 ,-iA_ 2>2 A 1>2 A_ 1>1 - Ai,_iA_i,ii4| i2 + Ax,-iA-x,i A -w + 

- 2 A_! -1^2,1^1,2^2,2 + 2 -42,-l-4 2 ,2-4i,2-4_i,i ~ 2 A X ,-1^2,1^-1,2^-2,2 - 

- 2 ^2,-1^2,2^1,1-4-1,2 + 2 A a ,-iA_a )2 -Ai,iA_i ) 2 + A\ 2 + 
+ 2 ^2,-i-4i,2vl_i,i + -4i,_i-4 2 ,i-4_i,2 - A_i,_iA 2 ,iAi,2 - 

- Al ) _lA_i i iA_2,2 - -4-2,2 + -4-l,-l-4i,iA_ 2 ,2 - A 2 ,-lAi,iA_i, 2 ^ cosfc 2 + 
+ (A^xA% 2 - -4-1,1 + 2 A2,lA_i, 2 A-2,2 - 2 -42,1-4-1,2-42,2 - 



Next, the coefficients in Eq.(|45|) defining asymptotical behavior of the generating functions of 
the return probabilities are given explicitly by: 



A [ ^_ x (u) = u 2 {u + 1) 2 ((tt - 2)u 2 - 2(tt 2 - 3tt - 2)u + n- 2), 
^i 2 -i( u ) = -ttu 2 (u + 1) 2 (u 2 + 2u + 2tt-3) x 

x ((2 7T - 3)u 2 + 2u+ 1)((tt - 2)u 2 + 4u + n - 2) 2 , 



A^l^u) = -7r(?i + l) 2 ((27r-3)M 2 + 2?i+l) 2 ((7r-2)ii 2 + 4M + 7r-2) 2 , 
A^^u) = u(7r-2)(u + l) 2 ((27r-3)u 2 + 2u+l), 

A^i^u) = -iru(u+l) 2 ((2Tr -3)u 2 + 2u+l)((ir -2)u 2 + 4u + tt -2) x 





+ (8 7T 2 - 34 7T + 52)lt 2 + (2 7T 2 - 8 7T + 16)ll + 7T — 2, 



x ((tt 2 - 4tt + 6)w 4 + (2 7T 2 - 6 7T + 8)u 3 - (2 7T 2 - 
- 207r + 28)u 2 + (27r 2 -67r + 8)u + 7r 2 -47r + 6) 
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A^l(u) = u 2 {{tt -2)u 4 + {2tt 2 -8tt + 16)u 3 + {8tt 2 -3Att + 52)u 2 - 
- (6tt 2 - 56tt + 80)u + 4tt 2 - 15tt+ 14), 

-U 2 TT (u + 1) 2 (U 2 + 2 U + 2 TT - 3) 2 ((7T - 2)M 2 + Au + TT - 2) 2 , 

(u + 1) 2 ((tt - 2)u 2 - 2(tt 2 - 3tt - 2) u + tt - 2), 
-7r(u + l) 2 (M 2 + 2it + 27r-3) x 
((2tt - 3)w 2 + 2u + 1)((tt - 2)u 2 + 4« + tt - 2) 2 , 
(tt - 2)u(?i+ 1) 2 (m 2 + 2u + 2tt - 3), 

-7r u(u + l) 2 (u 2 + 2 u + 2 tt - 3) ((tt - 2)u 2 + 4 u + tt - 2) x 
((tt 2 - 4 tt + 6)u 4 + (2 tt 2 - 6tt + 8) u 3 - 
(2 tt 2 - 20 TT + 28) u 2 + (2 tt 2 - 6 TT + 8) u + tt 2 - 4 TT + 6), 
u 2 (7r-2)(u+l) 2 (it 2 + 2M + 27r-3), 

-ttu 2 (u + l) 2 (u 2 + 2 u + 2 tt - 3)((tt - 2)?i 2 + 4u + tt-2)x 
{{tt 2 - 4 TT + 6) ?i 4 + (2 tt 2 - 6 TT + 8) u 3 - 
(2tt 2 - 20tt + 28) it 2 + (2tt 2 - 6tt + 8)u + tt 2 - 4tt + 6), 
(tt - 2){u + 1) 2 ((2 tt - 3)u 2 + 2u + 1), 

1) 2 ((2tt - 3)u 2 + 2 u + l)((?r - 2)u 2 + 4 U + 7T — 2) X 
((tt 2 - Att + 6 )u 4 + (2 tt 2 - 6 tt + 8) u 3 - (2 tt 2 - 20 tt + 28) u 2 + 

(2 7T 2 - 6 TT + 8) U + TT 2 - ATT + 6), 

-u{u 2 + (2tt - 2)u + l)~\u + l) 2 x 

((tt - 6)u 4 - 8 u 3 + (4 tt 3 - 16 tt 2 - 2 tt + 28) u 2 - 8 u + tt - 6), 
-ttu(u+ 1) 2 ((tt 2 -4tt + 6)m 4 + 

(2 TT 2 — 6 7T + 8)u 3 - (2tt 2 - 20?r + 28) u 2 + (2tt 2 - 6tt + 8)u + tt 2 - 4tt + 6) 2 , 
u{u 2 + (2tt - 2)u + 1)" 1 ((2tt 2 - 9 tt + 14) u 6 + 
(4 tt 3 - 20 tt 2 + 46 tt - 28) u 5 + (12 tt 3 - 66 tt 2 + 169 tt - 142) u 4 - 
(16 tt 3 - 168 tt 2 + 412 tt - 312) u 3 + (12 tt 3 - 66 tt 2 + 169 tt - 142) u 2 + 
(4 tt 3 - 20 tt 2 + 46 tt - 28) u + 2 tt 2 - 9 tt + 14), 

-TT U(U + 1) 2 ((7T 2 - 4 7T + 6) U 4 + (2 TT 2 - 6 TT + 8) U 3 - 

(2 tt 2 - 20 TT + 28) u 2 + (2 tt 2 - Qtt + 8) u + tt 2 - A TT + 6) 2 



X 

+ 



+ 



+ 



4%) 
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